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We study quantum-classical escape-rate transitions for uniaxial and biaxial models with finite spins 
S = 10 (such as Mni2Ac and Fes) and S = 100 by a direct numerical approach. At second- 
order transitions the level making a dominant contribution into thermally assisted tunneling changes 
gradually with temperature whereas at first-order transitions a group of levels is skipped. For finite 
spins, the quasiclassical boundaries between first- and second-order transitions are shifted, favoring a 
second-order transition: For Feg in zero field the transition should be first order according to a theory 
with S — > oo, but we show that there are no skipped levels at the transition. Applying a field along 
the hard axis in Fes makes transition the strongest first order. For the same model with S = 100 
we confirmed the existence of a region where a second-order transition is followed by a first-order 
transition [X. Martmes Hidalgo and E. M. Chudnovsky, J. Phys.: Condensed Matter (in press)]. 

PACS numbers: 75.45.+j, 75.50.Tt 



I. INTRODUCTION 

Recent experimental discovery of spin tunneling in 
large-spin compounds such as Mn^Ac and Fes (S = 10) 
also stimulated theoretical investigation of spin models 
with S ^> 1 which show different kinds of transition be- 
tween the classical mechanism of thermal activation over 
the potential barrier AU at temperatures To < T <§C AU 
and the quantum regimes involving tunneling under the 
barrier at T < Tq. The quantum-classical transition tem- 
perature To becomes well defined in the quasiclassical 
limit S — > oo and is of order To ~ AU/S, where AU 
is the barrier height. Possible types of the quantum- 
classical transition for a general model have been classi- 
fied in Ref. [j], the two main scenarios being the so-called 
first-order transition and the second-order transition. At 
the second-order transition the energy E* with which 
the system is crossing the barrier begins to move down 
from the value E c corresponding to the top of the bar- 
rier, and E*(T) approaches the bottom of the well E m [ n 
at T = 0. That is, for T < Tq there is a thermally as- 
sisted tunneling: Thermal activation up to the energy 
E = E* is followed by the tunneling at this energy level. 
At the first-order transition E* abruptly changes from E c 
to some lower value and then, again, E*(T) approaches 
the bottom of the well at T = 0. In this situation some 
interval of energy is skipped and it does not contibute 
to the escape from the metastable well at any tempera- 
ture. There are more exotic cases such as a second-order 
transition followed by a first-order transition. 

Whereas for quantum particles in the well it is diffi- 



cult to realize transitions other than the second order, it 
has been recently shown that there are all the types of 
quantum-classical transitions in the spin model with the 
Hamiltonian 



n 



-DS, + BS X — H X S X 



HySy 



H Z S Z 



(1) 



which is convenient to parametrize in terms of the re- 
duced hard-axis anisotropy b = B/D and the reduced 
fields h x = H X /(2SD), etc. In Ref. | it was shown that 
in the uniaxial model (B = 0) transverse field controls the 
order of transition which is first order for small transverse 
fields. In Ref. || the exact quasiclassical value of the crit- 
ical transverse field h — h c — 1/4 has been obtained. In 
Ref. ^ the whole phase diagram of escape-rate transitions 
for the uniaxial model in the plane h x , h z was drawn, the 
boundary line h xc (h z ) going to zero at h z = 1. 

For the biaxial model, B =/= 0, in zero field the transi- 
tion was shown to be first order for b < b c — 1 and second 
order for b > 1.0 Longitudinal field suppresses the value 
of b c , so that b c (h-) vanishes at h z = 1 together with the 
potential barrierfl The quasiciaasical result in this case 
reads b c = (1 - h 2 z )/(l + 2h 2 z )W 

Field along the medium axis, H y , in Eq. [j] also favors 
the second-order transition. Some points of the boundary 
between the first- and second-order transitions have been 
obtained in Ref. || numerically, whereas the analytical 
expression for the boundary was obtained later in Ref. ||. 
Recently the phase diagram of the biaxial model with the 
fields along medium (H v ) a nd easy (H z ) directions has 
been considered in Ref. [l(j. The qualitative results are 
a combination of those of Ref. ^ for the uniaxial model 
and those of Refs. 0, || for the biaxial model with the field 
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H z : Increasing of all b, h y , and h z favors the second-order 
transition. 

The most interesting model is the biaxial model with 
the field along the hard direction H x , in which oscilla- 
tions of the tunneling probability as function of h x haiza 
been established theoreticallyEil and experimentally.^ 
The phase diagram of escape-rate transitions for this 
model has been recently considered in Refs. |l3|, Un- 
like all other models, there is a first-order transition for 
h x ~ 6, even for the large enough values of h x and b which 
would alone cause a second-order transition. That is, the 
order of transition can change as II-I-II with increasing 
h x or b. Moreover, in Ref. [l4] ranges of parameters have 
been found where the second-order transition is followed 
by a first-order transition or a first-order transition is 
followed by another first-order transition with lowering 
temperature. 

Thus the theoretical investigation of the escape-rate 
transitions in the spin system described by Eq. [I] in 
the quasiclassical approximation 5 3> 1 is nearly com- 
pleted. Experimentally studied materials, however, have 
the moderate spin value 5 = 10, which can result in de- 
viations from the predictions of the quasiclassical theory. 
Indeed, in Mn^Ac in zero applied field one can expect a 
stpang first-order transition but experiments of Kent et 
atB show only one skipped energy level, m = —9. For 
Fes (b = 0.47) in zero field one expects a first-order tran- 
sition but recent measurements of WernsdorferE3 suggest 
that each energy level becomes dominant in the escape 
at some temperature, i.e., the transition is second order. 

Since it is very difficult to find 1/5 corrections to the 
quasiclassical results, one has to look for alternative ap- 
proaches. For finite spins, the problem can be solved 
in a purely numerical way, and the calculations can be 
performed on a modern PC within a reasonable time for 
5 < 100. This is the aim of the present article — to find 
out which energy levels make the dominant contribution 
to the escape rate at different temperatures for different 
particular cases of the spin model with the Hamiltonian 
of Eq. [|. We show that, in accordance with predictions 
of quasiclassical model, the dominant level docs not nec- 
essarily change continuously from the top to the bottom 
of the mctastable well when temperature is lowered. 

The rest of the article is organized as follows. In Sec. |l| 
we reformulate the theory of thermally assisted tunneling 
in terms of quantities which can be directly computed for 
finite spins. In Sec. ?? we consider the uniaxial model 
with transverse and longitudinal fields and make a com- 
parison of exact numerical results for the temperature 
dependence of the tunneling level with earlier perturba- 
tive results. In Sec. IV the calculations are performed for 
the biaxial model with the field along the hard direction. 
We confirm the existence of more complicated scenarios 
of the escape-rate transitions for this model. 



II. BASIC FORMALISM 

Quasiclassical approach to the quantum statistical 
metastability considers the spectrum of quantum states 



as continuous and uses the following expression for the 
escape rate (see, e.g., Refs. [It], 



r~ JdEW{E)e 



(B-B min )/T 



(2) 



where W{E) is the probability of tunneling at an energy 
E. The latter can be written as 



W(E) 



1 



1 + exp[S{E)] ' 



(3) 



where for the barriers parabolic near the top the 
imaginary-time action S(E) goes linearly through zero 
for E crossing the barrier top level E = E c and it is 
analytically continued into the energy region above the 
barrier. In the latter case formula (|3|) describes quan- 
tum reflections for a particle going over the barrier, with 
W(E) slightly lower than 1, whereas for the energies be- 
low the top of the barrier W(E) is exponentially small 
in the quasiclassical case. The action S(E) can be calcu- 
lated for spin systems with a number|-af-different meth- 
ods such as the instanton approach Jl§HP | 'jl mapping on a 
particle with the WKB apprGecimation^'EfH'Bi and the 
discrete spin WKB method.Ea 

At higher temperatures the integral in Eq. @ is domi- 
nated by E ~ E c which results in the Arrhcnius temper- 
ature dependence of the escape rate T = Tq exp(AC//T), 
where AC/ = E c — E m - m . At lower temperatures the rel- 
evant region of energies goes down, which is the regime 
of thermally assisted tunneling. Since for quasiclassical 
systems the crossover between the two regimes occurs 
at a temperature T -C AC/, the integrand in Eq. (||) is 
a product of two very rapidly increasing and decreasing 
functions of energy and thus can be approximated by 



T - max;; W{E)e 



-(E-E miB )/T 



(4) 



Within this approximation, the crossover between the 
classical and quantum regimes becomes a transition at a 
well defined temperature To. The mathematical descrip- 
tion of this transition is analogous to the well known phe- 
nomenological Landau model of phase transitions (the 
Landau theory), as was pointed out in Ref. ^. The tran- 
sition can be second or first order. It should be stressed, 
however, that it is only a formal analogy and there are 
certainly no many-body effects in the problem of the es- 
cape rates we are studying. Integration across the maxi- 
mum of the integrand in Eq. (Jfy smears the transition and 
transformes it into a crossover. In the case of a second- 
order transition the width of the crossover region around 
To is AT cx 1/y/S (Ref. ^) and disappears in the quasi- 
classical limit. In the case of a first-order transition there 
are two competing maxima of the integrand in Eq. (Q), 
and the transition is from one maximum to the other. In 
this case the width of the crossover region is even smaller: 
AT cx 1/5 Jj We should stress that in spite of the smear- 
ing of the escape-rate transition, there is always a funda- 
mental difference between the two situations: one shift- 
ing with temperature maximum of W(E) exp(— E/T) or 
two competing maxima of W(E) exp(—E/T). We will 
illustrate this difference for various models below. 
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It is convenient to express the tunneling probability 
through the quantities which can be directly computed, 
using the quasiclassical formula for the tunnel splitting^ 



7T 



S(E) 



(5) 



where loe is the frequency of the oscillation in the well 
at the energy E. Using this formula one obtains 



W{E) 



(AE)' 



(6) 



{w E /K) 2 + {AEy 

Since, again, AE given by Eq. (|J) becomes formally much 
larger than uje for the energies above the top of the bar- 
rier if S(E) is analytically continued into that region, Eq. 
(H) gives W(E) fast approaching 1. 

Eq. (^) is the starting point for the numerical solu- 
tion of the problem for finite spins. Wc consider the 
situation where there are pairs of quasi-degenerate lev- 
els in different potential wells and we compute the tun- 
nel splittings for these pairs numerically. The oscilla- 
tion frequency loe is nothing else than the difference of 
the energy of the adjaicent levels in one of the wells: 
uAE n ) = 6E n = E n+ i — E n . The discrete analog of Eq. 
(|is 



(AE n ) 2 6E n 



^ (6E n /ir) 2 + (AE n ) 2 ' 



-(B„-E min )/T 



For large spins and low temperatures one has 



r~max„exp[-F(£;„)/T], 



(J) 



(8) 



where 



exp 



F{E n ) 



T 



(AE n ) 2 6E n 



[5E n /n) 2 + (AE n y 



-(-E„--E mi „)/T 



(9) 



We will call the energy level minimizing the effective free 
energy F (Ref. |^) and thus maximizing the combined 
probability of escape the tunneling level or the level of 
thermally assisted tunneling. 

For the energy levels below the top of the barrier, one 
has AE <C SE n , whereas above the top of the barrier 
the levels are not grouped in pairs, i.e., formally, AE ~ 
5E n . Since for large spins the transition between the two 
ranges of energy is rather sharp and below the top of the 
barrier AE changes much faster-than 5E, one can look 
for the maximum of the functiona 



(AE T 



)2 e -(E„ 



-Emin)/T 



(10) 



Although above the barrier one cannot strictly speak of 
tunneling, the formula above gives correct results since 
AE n becomes weakly dependent on energy and this re- 
gion is suppressed by the fast decreasing Boltzmann ex- 
ponential. In this paper, we will use Eq. (|9|) instead of Eq. 
( |l0| ) since we are going to make a comparison between 
the exact mumerical solution and the solution that uses 
the perturbative formula for the level spittings .cj Since 
the latter gives AE ^S> SE n above the barrier, Eq. (||) 
is more appropriate because it gives physically correct 
results in this energy range. 



III. UNIAXIAL MODEL WITH EXTERNAL 
FIELD 

This model is the first of spin models for which 
the first-order escape-rate transition has been found 
theoreticallyn in the region of small transverse fields 
H x using .. tijie perturbative formula for the level 
splittingsr 4 ' ~ 



Ae, 



2D 



[(to' 



1) 



112 



' [S + m')\{S - m)\ (H*_ 
(S-mI)\{S + m)\ \2D 



(11) 



Here the longitudinal field enters through the resonance 
condition 

E rn = E m > , to < 0, to' = —to — k 

H z = H zk = kD, fc = 0,±l,±2,... (12) 

with E m — —Dm 2 — H z m, see Fig. 1 of Ref. [|. Recently 
corrections to this formula have been ontained in Ref. 
p6| . The use of the perturbative formula for the splittings 
cannot, however, give an acccurate value of the boundary 
h xc = H xc / (2SD) between the first- and second-order 
transitions since the transition occurs at h x — h xc = 1/4 
(Rcf. ||) which is not small. 

Our next task is to perform a purely numerical calcu- 
lation illustrating first- and second-order transitions for 
finite spins in the transverse field of arbitrary strength. 
For H x ^ 0, spin projections on the z axis, to, are 
no longer good quantum numbers. We will continue, 
however, to enumerate the exact levels in terms of to 
to keep a link to the previous work. At the /cth res- 
onance, the lowest k levels are not splitted and local- 
ized in the right well. We will formally ascribe them 
S z values to = S, S — 1, . . . , S — k + 1. Higher levels 
are grouped in tunnel-splitted pairs which we denote as 
{to, to'} = {-S, S - k}, {-S + l,S-k- 1}, etc. 

For the diagonalization of the spin Hamiltonian we 
used Wolfram Mathematica which allows one to perform 
calculations with any desired precision. We used the pa- 
rameter set of Mni2 Ac and ignored the anisotropy of the 
type D4S* for simplicity. Our numerical results for the 
level splittings reproduce those of Ref. [27|, where a quan- 
tum dimcr problem, which is mathematically identical to 
the spin-in-field problem, has been studied. 

The perturbative and exact results for the temperature 
dependence of the tunneling level mx which maximizes 
Eq. (S) are shown in Fig. [I] for H z = 0. For the field 
h x = 0.125 which can be considered as small, the pertur- 
bation theory well describes the transition temperature 
To and the order of transition. The only noticeable dis- 
agreement with the exact results is that regarding the 
hight of the barrier which is visualized here through the 
value of tot in the classical regime. This is not a surprize 
since the PT breaks_dQwn near the top of the barrier for 
whatever small /i^E^El For h x = 0.125 and S = 100, 
many levels are skipped at the transition temperature, 
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1. Perturbative and exact results for the temperature 
dependence of the tunneling level mi for the uniaxial spin 
model with transverse field. One can see that the PT holds 
for small h x and that finite values of the spin 5* favor the 
second-order transition. 



thus this transition is first order. For h x = 0.125 and 
5 = 10, the skipped range is smaller, and the situation is 
closer to a second-order transition than that for S = 100. 
The dependence on the spin value is even more clearly 
seen for h x = 0.25 which is the exact boundary between, 
first- and second-order transitions in the limit S — * oo.u 
For S = 10 there are no skipped levels and the depen- 
dence rriT(T) is far from a jump. For S — 100 there are 
no skipped levels, too, but the dependence m<r(T) has a 
rather high slope near To. In the limit S — > oo, the low- 
slope part of the dependence itlt(T) at T > Tq becomes 
horizontal, and the derivative dmr/clT becomes infinite 
at T = To — 0. For S = 100, there are no skipped levels 
even for h x = 0.2. 

On Fig. H we show exact numerical results for mr(T) 
for the S = 10 and S = 100 models with h x = 0.125 and 
two values of the longitudinal field, h z = and h z = 0.4. 
These results confirm that increasing of h z drives the 
system into the region of the second-order transitions.u 
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FIG. 2. Temperature dependence of the tunneling level rriT 
for the uniaxial spin model with transverse and longitudinal 
field. 
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FIG. 3. For the model H = -DS 2 Z + BS 2 X - H X S X phase 
diagram includes regions of first- (I) and second-order (II) 
escape-rate transitions, as well as the regions where a sec- 
ond-order transition is followed by a first-OKl^r one (II-I) or 
the regions of the transitions of the I-I type.l 



IV. BIAXIAL MODEL 



We will concentrate here on the mos 
with the field along the hard directionl 



BBS 



ng model 



H = -DSl + BS 2 ~H x S x . 



(13) 



The phase diagram for the model above, which has been 
obtained in Ref. |lj is shown in Fig. ||. The boundaries 
marked a and b have also been obtained in Ref. 13. Apart 
from regions of the first- and second-order transitions 
marked by I and II, this phase diagram contains the re- 
gion where a second-order transition is followed by the 
first-order one (II-I) and a rather narrow region where 
a first-order transition is followed by another first-order 
transition (I-I). The possibility of such multiple transi- 
tions has been predicted in Ref. [I] and here is their first 
realization in a spin model. 
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FIG. 4. Temperature dependence of the tunneling level rrir 
for the biaxial spin model with b = B/D = 0.470 in zero field. 
Note the second-order transition for S = 10. 



Let us now draw the plots of my(T) for S = 10 and 
S = 100 for different transverse fields h x for the value of 
the transverse anisotropy b — 0.47 which is appropriate 
for Fes- In zero field one expects a first-order transi- 
tion for b < 1 in the quasiclassical limit B This is con- 
firmed by our results for S = 100 in Fig. |. However, 
for S = 10 there are no skipped levels, although itit(T) 
goes rather steep. This is in accord with recent experi- 
ments by Wernsdorfer on -Ees in zero field, which suggest 
a second-order transitionJl3 

The behavior changes strikingly if a sufficiently strong 
field h x is applied. One can see from the Fig. || that for 
h x = 0.5 for both S = 100 and S = 10 the transition is 
the strongest first order. This effect should be observable 
on Fes- Further increasing the field makes the potential 
wells so shallow that there are only few levels left. This 
makes difficult to make a judgement about the order of 
the transition for S = 10. For h x = 0.75 there are no 
skipped levels for the S = 10 model and one could speak 
about a second-order transition. For S — 100 one can 
clearly see a second-order transition followed by a first- 
order transition, in accordance with the phase diagram 
on Fig. |[ Dependences rriT(T) for S = 100 and many 
different values of h x are shown in Fig. |^. The different 
types of transition in Fig. ^ are in accord with the phase 
diagram of Fig. a. 

The most exotic behavior of itit(T) takes place for 
small values of transverse anisotropy and field, where one 
expects two first-order transitions (see Fig. ||). The be- 
havior of m T (T) for S = 100, b = O.OlJ^and h x = 0.1 in 
Fig. ^ confirms the prediction of Ref. |lj and shows two 
jumps. Note that with lowering temperature rriT{T) for 
S = 100 begins to go down continuously and then makes 
the first jump. Thus one could speak about the succes- 
sion of transitions of the type II-I-I, where the second- 
order transition is solely due to the finite value of the 
spin and vanishes in the quasiclassical limit. The same 
effect also takes place in a simpler uniaxial model with a 
transverse field h x < h xc . For 5 = 10 the jump at higher 
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FIG. 5. Temperature dependence of the tunneling level mr 
for the biaxial spin model with b = 0.470 in hard-axis fields 
h x = 0.5 and 0.75. 



temprerature disappears and one thus has a second-order 
transition followed by a first-order transition. We illus- 
trate the behavior of the effective free energy F of Eq. (||) 
for b = 0.015 and h x — 0.1 in Fig. |[ (For convenience, 
we use the value SD = 2.34 K of Fes.) Since the de- 
pendence F(T) on the energy level is extremely flat near 
T = 0.21 K, the approach using Eq. (||) instead of Eq. 
(H) is valid for rather high values of S. 

As we have mentioned in the Introduction, for the bi- 
axial model with the field along the hard axis and the 
integer spin, tunneling is quenched wheneverliU 



H x = (1 + 2n)y/B(B + D), 
n = -S,-S + 2,...,S-l. 



(14) 



In the quasiclassical formalism it manifests itself in the 
vanishing of the prefactor in the tunneling probability. 
For S ^> 1 the role of the prefactor is difficult to see 
because the exponential terms dominate. The rate of 
thermally assisted tunneling (in the log scale) and thus 
the transition temperature To, which depends logarith- 
mically on the prefactor, are significantly reduced only in 
very close vicinities of quenching points. For moderate 
spins such as S = 10, the quenching effect may be quite 
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FIG. 6. Temperature dependence of the tunneling level mr 
for the biaxial spin model with b — 0.470 and S = 100 in 
different hard-axis fields. 
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FIG. 7. Temperature dependence of the tunneling level hit 
for the biaxial spin model with b = 0.015 and h x =0.1 show- 
ing two first-order transitions for S = 100. 
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FIG. 9. The escape rate F in Fes vs 1/T for different trans- 
verse fields. For h x — 0.475 tunneling is almost quenched. 



substantial. On Fig. |9] the tunneling rate and the value 
of To are suppressed near h x = 0.457 which corresponds 
to n = 5 in Eq. (Q). 

V. DISCUSSION 

Our direct numerical investigations of quantum- 
classical escape-rate transitions in spin models with fi- 
nite S confirmed predictions of quasiclassical approaches 
in the case of large S and revealed deviations to the favor 
of a second-order transition for moderate spins. In partic- 
ular, in Feg in zero field the moderate spin value S = 10 
makes the transition second order. On the other hand, 
applying a field along the hard anisotropy axis makes 
the transition in Fes the strongest first order, which can 
be probably observed in experiment. For some values of 
the field tunneling is quenched and the rate of thermally 
assisted tunneling drops down. 



For the biaxial model with the field along the hard 
axis, we numerically confirmed the existence of the re- 
gions where (i) a second-order transition is followed by 
a first-order transition and (ii) a first-order transition is 
followed by another first-order transition with lowering 
temperatureJU This model seems to be the only model 
up to date which demonstrates such a complicated be- 
havior. 
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